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Abstract 

In recent work of Baik, Deift and Rains convergence of moments was established 
for the limiting joint distribution of the lengths of the first k rows in random Young 
tableaux. The main difficulty was obtaining a good estimate for the tail of the distribu- 
tion and this was accomplished through a highly nontrival Riemann-Hilbert analysis. 
Here we give a simpler derivation. The same method is used to establish convergence 
of moments for a random growth model. 

Since the paper of Baik, Deift and Johansson [0J in which the authors determined the 
limiting distribution of the length of the longest increasing subsequence in a random permu- 
tation, or equivalently the length Ai of the first row in a random Young tableau, there have 
been a variety of extensions and generalizations. In [§] the same authors determined the lim- 
iting distribution of the length A2 of the second row and Borodin, Okounkov and Olshansky 
in || and Johansson in [[TIJ determined the limit of the joint distribution of Aj, • • • , for 



any k. The result was that for fixed xx, • • • , x k the limit 



lim Pr (\i < 2VN + x f Y 1/6 , i = 1. 



1 k 

N 



exists and equals the limiting joint distribution of the largest k eigenvalues in the Gaussian 
unitary ensemble of random matrices. (Here N was the number of boxes in the Young 
tableaux, which were given the Plancherel measure.) 

In the end it turned out to be the question of the asymptotics of certain Fredholm 
determinants, and their derivatives, associated with the operator K n acting on £ 2 (n, n+1, • • •) 
whose matrix entries are given by 

K M = ±fc) (g) . (!) 

k=l K^+Ji+k \V-J-k-j 

Subscripts denote Fourier coefficients and here (p±(z) = e tz±1 , the Wiener- Hopf factors of 

<p(z) = e****' 1 ). 

In order to prove convergence of moments of the distributions one has to have uniform 
estimates for the probabilities and the only difficulties occur at —00. Convergence of moments 
was established for Ai and A2 in [I], |2| and, for the joint distribution, in [|J. Essentially what 
one has to show is that each 

d^ 

— det(I-vK n ) 
av v=l 
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is very small if 2t — n is large. In the last cited paper an integral operator was introduced 
which was an alternative to K n and which was amenable to a highly nontrivial Riemann- 
Hilbert analysis which gave the necessary estimates. Here we shall show that a general 
determinant inequality allows one to derive the estimates more easily. 

In H the authors studied a random growth model called oriented digital boiling and 
determined the limiting shape and its fluctuations. In this case the probability in question, 
Pr (H < h) where if is a certain random variable, is equal to det(J — K^) where now 

V (z) = (l + zT(l-rz- 1 )- m , 

with Wiener-Hopf factors 

¥> + (z) = (1 + ^(z) = (l-rz- 1 )- m - (2) 

The asymptotics were determined in various regimes. In the GUE universal regime of M 
n = am with a fixed and ar < 1. It was shown that for a certain constant c, the "time 
constant", the limit 

lim Pr (H < cm + s m 1 ^ 3 ) 

m— >oo 

exists and equals the limiting distribution function for the largest eigenvalue in the Gaussian 
unitary ensemble. To establish convergence of moments the main problem is to show that 
the determinant is very small, uniformly in m, for large negative s. (Estimates established 
in [0] show that the probability approaches 1 exponentially for large positive s, uniformly in 
m.) 

The inequality which allows us to obtain adequate estimates easily is the following. 
Lemma 1. If K is a trace class operator all of whose eigenvalues are real and < 1 then 

det{I-K) < e~ tlK . 



Proof. If A < 1 then < 1 - A < e~ A . Take the product over all eigenvalues of K. 

We first apply the lemma to the Young tableaux distributions. As observed in Q, since 
all eigenvaues of K n in this case are real and lie between and 1, it is enough to estimate 
det(l — vK n ) for any fixed v G (0, 1). If 2t — n + sn 1 ^ 3 the authors obtained bounds of the 
form 

e~ Ss3/2 if s<n 2 / 3 



det(I -vK n ) < 



e- st if s > n 2 ' 3 



for some 5 > 0. (These are equivalent to inequalities (5.10) and (5.11) of M.) These bounds 
will follow from the corresponding estimate for the trace. 

Proposition 1. We have for some 5 > 

5s 3 / 2 iis<n 2 / 3 



trK n > 



St if s > n 2/3 . 
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Proof. We have 

k=l \f+J n+ k \<P+J - n -k 

and here (ip^/ip + )k = (<p+/<p-)-k — ^(2t), where Jk is the bessel function. Therefore 

oo 

tiK n = J2k\J n+k (2t)\ 2 . (3) 
k=l 

The uniform asymptotics for the Bessel function (J7|, sec. 4.8) give for u > 1 

j fc (A;u) = Qjfc 2 /W(tt)) V2 Ai(-A; 2 / 3 ^(u)) (1 + Oik' 1 )), 

where ^ip(u) 3 ^ 2 = /"(l — t -2 ) 1 / 2 (it. (The error term has to be modified near the zeros of the 
Airy function, that is, when ~kip(u) 3 ' 2 is near the set — | + 7rZ.) It is easy to see that if 

2t = k + s k 1/3 

then k 2 l 3 ip{2t/k) tends to +oo as k, s — > +oo. Hence, from the asymptotics of the Airy 
function we have for large k and s and some constant 5 > 

wn>* r*&(%)Y% ftY" (4) 

as long as ^kip{2t/k) 3 ^ 2 is bounded away from — j + nZ. 
Now t and n are related by 

2t = n + sn 1/3 . 
We replace k by n + k in (§) and choose those such that 

«isn 1/3 < jfe < a 2 sn 1/3 , (5) 

where ai < ol^ < 1. Then we find that 

, . sn~ 2 ^ 3 if s < r?l 3 

- 1 



n + k [l iis>n 2 / 3 , 

where "~" means "is bounded above and below by positive constants times" . From this and 
(|) we deduce that with another 5 > 



f 5s- 3 / 4 if s < n 2 ' 3 



\Jn+k{2t)\ > { 



Sn - 1/es -l/2 i{s>n 2/3 



2 We should be careful because of the requirement that |(n + k)ip(2t/(n + fc)) 3//2 be bounded away from 
— j + 7rZ. If a.\ in (0) is very close to 1 then two consecutive values of (n + k)ip(2t/(n + fc)) 3 / 2 differ by at 
most a small constant. Thus with such a choice of ot\ we can be assured that the estimates hold for a large 
fraction of the k satisfying (||). 
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And from this, summing in (^) over only those k satisfying (BJ), we deduce that for s < n 2 ^ 3 
the trace of K n is at least a constant times s 1//2 n 2 / 3 which is in turn at least a constant times 
s 3 / 2 and for s > n 2 ^ 3 it is at least a constant times sn 1 ^ 3 which is in turn at least a constant 
times t. This proves the proposition. 

We turn now to the growth model where Kh is defined in terms of the Wiener-Hopf 
factors (D-0 

Lemma 2. The eigenvalues of K% are all real and less than 1. 
Proof. An expression for K(i,j) when 

cp(z) = (i+tzr(i+tz- i )- m 

is given by formula (4.9) of Pfl The i,j entry of our operator is equal (zr 1 / 2 )* - - 7 times the 
i,j entry of this one in which £ = ir 1 / 2 . From this and the formula given there we see that 
our K(i,j) is of the form a* L(i, j) bj where a;, bj > and L(i,j) is symmetric. Hence all 
eigenvalues of Kh are real. Next, notice that det (J — Kh) is always positive, being a nonzero 
probability. So no eigenvalue of Kh can equal 1. The eigenvalues being continuous functions 
of r, and the eigenvalues being all when r = 0, we see that no eigenvalue can be > 1. 

Recall that n = am with a fixed, ar < 1 and 

Pr (H < cm + s m 1/3 ) = det (J - K h ) 
where h = cm + sm 1 / 3 . The time constant c is given by the formula 

c = — !— ((1 - a)r + 2Var). 
1 + r 

Since the only problem occurs when s is large negative we change s to — s and think of s as 
large and positive. So we have to estimate det (/ — Kh) when h = cm — s m x l 3 . By Lemmas 1 
and 2 it suffices to find a lower bound for the trace. 

Lemma 3. If s is sufficiently large and sm~ 2 l 3 sufficiently small there exists a constant 
5 > such that tr K h > 5s 3/2 . 

Proof. Given functions / and g we have 

X! fi+k 9-k-j = —T~2 I I f( Z ^ 9{Z2) X Z l l ~ k ' 1 Z 2 +3 ^ dz l dz 2 
k=l k=l 

3 Now the operator is not trace class on the usual £ 2 space. From the integral formula for the matrix 
entries given at the beginning of the proof of Lemma 3 one sees that they grow in i as i n_1 and decrease 
expontially in j as j — * oo. So we must take our £ 2 space to be weighted. If the weight function is w(i) = 2~ 4 , 
for example, the operators Kh are trace class. 

4 Our K(i,j) differs from that of |j| by a factor (— This does not affect the eigenvalues and so we 
make believe the factor is not there. 
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where the integrals are taken over circles about the origin on which \z 2 \ < \z±\. If we 
consider these as matrix elements of an operator acting on £ 2 (h, h + 1, ■ ■ •) then the trace of 
the operator equals 



-h Ji 



^ 2J j fWgtojj—^dKdz,. 
In our case we set 

f(z) = (1 + z)- n (1 - rz- 1 )-™, g{z) = (1 + z) n (1 - rz~ x ) m 
and define, following the notation of ||, 

iP(z, d) = (1 + z) n (1 - rz- v ) m z c ' m . 



If h = dm we have 



Iff ip(z 2 ,c') dz l dz 2 

tr K h ^ 



47r 2 J J ip(zi,d) (zi - z 2 ) 2 ' 

Here we should have \zi\ < 1 and \z 2 \ > r, so r < 1. The formulas for the matrix entries and 
trace hold for all r by analytic continuation if the z 2 contour is inside the Z\ contour and the 
latter has r on the inside and —1 on the outside; the z 2 contour has on the inside. 

Now h = cm — sm 1 / 3 , so d = c — sm^ 2 ^ 3 , and we want a lower bound for the trace. We 
use steepest descent on each of the two integrals as if the factor (z\ — z 2 ) ~ 2 were not there. 
When d = c the function 

a(z, c) = — log ip(z, d) = a log (1 + z) + log (r — z) + (c — 1) log (—z) 
m 

has vanishing first and second derivatives at the single critical point u c G (—1, 0) and the two 
steepest descent curves are as shown in Fig. 6 of p. The one for the z 2 integration, denoted 
here by C 2 , is closed and contains on the inside while C%, the one for the Z\ integration, has 
two arms going to infinity and separates from —1. Our original contours can be deformed 
to these. 

For d c there are in general two critical points uf, given by 



-(1 - r)d - (1 - a)r ± J((l + r)d + (a - 1) r) 2 - Aar 

2 (a + c) 

The square root is zero when d = c. When d decreases from c the square roots are, at 
least for a while, purely imaginary^ Denote by u^, the one in the upper half-plane and by 
u~i the one in the lower half-plane. As we can see the steepest descent contours C\ and C 2 
now cross and so C 2 is no longer inside inside G\. (Also, the resulting double integral is not 
absolutely convergent but must be interpreted as a principal value). Therefore our original 

5 There are cases when the square roots are real for some c' < c, and the following descussion would have 
to be modified when this happens. See Remark 1 following the proof of the lemma. This is the reason for 
our assumption that smT 2 / 3, = c — c! is sufficiently small. 
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double integral is not equal to the integral over C\ xC%. In fact, we can see that our original 
double integral (without its factor — 1/4tt 2 ) is equal to the integral over G\ x C 2 plus 

/ r ip(z 2 ,d) dz 1 dz 2 
vT, Jv ip(z!,d) (zi - z 2 ) 2 ' 

where the z 2 path of integration is to the right of —1 and to the left of and the loop V 
enclose all point of the z 2 path. The inner integral is equal to 



ip(z 2 ,c') 

and integrating this and dividing by the factor — Att 2 gives 

I Tfl 

7 [log V> fee') - log -0 fee')] = -— [a(u^, d) - a(uj,c')], 



2ni l orv c ' ' 2irV 

where in the computation of the difference of logarithms we go over a path to the right of 
— 1 and to the left of 0. 

It follows from the displayed formula after (3.17) in H that 

re' 

c ') - c ') = / l°g -r rf 7- 

Jc 

The logarithm in the integral is positive purely imaginary and d < c, so the left side is 
negative purely imaginary. Using (|6]) we see that in fact this expression times its factor 
—m/2ni is at least a positive constant times m(c — c') 3 / 2 = s 3//2 . 

It remains to estimate 

ip(z 2 , d) dzi dz 2 



lc 2 Jci ij){zi,d) {z x - z 2 ) 2 ' 

and we shall show that it is 0(1). It can be seen that C\ closes at —1 while C 2 closes at 0. 
Denote by Cf the portions of these contours in the upper half-plane and by the portions 
in the lower half-plane. The distances between Cf and C 2 are of order (c — c') 1 / 2 . Lemma 6 
of H tells us that a"{u^,,d) is of the order ±(c — c') 1 / 2 . It follows that the integrals over 
Cf x C 2 are of order (c — d)~ 1 m~ 1 (c — c') -1 / 2 = s~ 3//2 . So these integrals are 0(1) for s > 1 
and we may restrict attention to the integrals over Cf x C 2 . 

Considering the integral over x C 2 we can confine attention to an arbitrarily small 
neighborhood of (u^,,u^,). (The integral over the complement of any such neighborhood is 
exponentially small.) For notational convenience we write u for uj. In our neighborhood we 
have 

C ') = e m[ CT "(«,c')((^-«) 2 -(^i-«) 2 )+0(| 22 -«| 3 +^i-«| 3 )] 

ip(zi,d) 

Because our contours cross at a positive angle we have (z{ —u) 3 /^ — z 2 ) 2 = 0(\zi —u\) and 
it follows that the contribution of the O term on the right is at most a constant times 



(mV2( c - c') 1 / 4 ) 3 
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So this contribution is also 0(1) for s > 1. 
That leaves the main integral 

gtiw" (u,c')((z2-u) 2 -(zi-u) 2 ) 



C+ (zi - Z 2 ) 2 



dz\ dz 2 



which, recall, is interpreted as a principal value. Recall also that we confine attention to a 
neighborhood of (u,u). If we make the variable changes z% — u = QJ (ma" (u, c')) 1 / 2 with an 
appropriate choice of square root the above becomes 

(C|-C?) 

d(i d( 2 , 



r 2 M (Ci - C: 

where Tj are long contours crossing at (0,0), with Ti horizontal there and T 2 vertical. The 
integral outside a neighborhood of (0, 0) is bounded. In a neighborhood of (0, 0) we can 
replace the exponential by 1 with error 0(1) and the principal value integral resulting after 
this replacement is 0(1). This establishes the lemma. 

Remark 1. We know that c is a value of d such that there is only one critical point. But 
if 2^far < (1 — a)r there will be another nonnegative value of d where there is only one 
critical point, namely 

-((1 — a)r — 2y / ar). 



1 + r 

And for smaller d there will be two real critical points. So the nature of the curves Cj will 
be different. In particular for d less than this value C 2 is outside C\. What happens in 
these cases is that instead of integrating from u~, to u^, in (|7|) the outer integral is a closed 
loop with the result that the integral with its factor — 1/Att 2 is equal to n exactly. However, 
concerning ourselves with this is unnecessary since we know enough already. 

Proposition 2. There exists a 5 > such that det(J — Kh) < e~ 5s3/2 for sufficiently large s. 

Proof. For sufficiently small sm" 2 ' 3 , say sm~ 2 ^ 3 < 77, Lemmas 1-3 give the stated estimate. 
In particular for sm~ 2 ^ 3 = rj the determinant is at most e~ Ss3/2 = e ~ Sr i 3/2rn , Since the 
determinant is a nonincreasing function of s this bound holds also for sm~ 2 ^ 3 > t]. But since 
h = cm — sm 1 ' 3 > we have m > c~ 3 / 2 s 3 / 2 and therefore the stated bound holds also for 
these s, with 5 replaced by S(r]/c) 3 ^ 2 . 

Remark 2. Baik, Deift, McLaughlin, Miller and Zhou have recently proved convergence 
of moments, using Riemann-Hilbert techniques, for a growth model of Johansson |J. There 
an anologous determinant arises where K(i,j) is given by (P with 

y? + (z) = (l+tz) M ^(z) = (l+tz-Y. 

We thank the authors for alerting us to their results. Very likely their methods could also 
be used for our growth model. 

Remark 3. The proof of Lemma 2 is not very satisfactory. There should be a reason for 
its truth and it should be a special case of a more general fact. There may be some merit in 
the following 



7 



Conjecture. Let r« and Sj be finitely many nonegative real numbers. Then the eigenvalues 
of the operator K n acting on £ 2 (n, n + 1, • • •) whose matrix entries are given by (D with 



z 



<p+(z) = n c 1 + r ^)' w = n (! - ^ 

* i 
are real and lie between and 1. 

Numerical evidence supports the truth of the conjecture. 
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